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Abstract: We construct models of dark matter with suppressed spin-independent scat¬ 
tering cross section utilizing the existing simplified model framework. Even simple com¬ 
binations of simplihed models can exhibit interference effects that cause the tree level 
contribution to the scattering cross section to vanish, thus demonstrating that direct de¬ 
tection limits on simplified models are not robust when embedded in a more complicated 
and realistic framework. In general for fermionic WIMP masses > 10 GeV direct detection 
limits on the spin-independent scattering cross section are much stronger than those coming 
from the LHC. However these model combinations, which we call less-simplihed models, 
represent situations where LHC searches become more competitive than direct detection 
experiments even for moderate dark matter mass. We show that a complementary use of 
several searches at the LHC can strongly constrain the direct detection blind spots by set¬ 
ting limits on the coupling constants and mediators’ mass. We derive the strongest limits 
for combinations of vector -I- scalar, vector -|- “squark”, and “squark” -|- scalar mediator, 
and present the corresponding projections for the LHC 14TeV for a number of searches: 
mono-jet, jets -|- missing energy, and searches for heavy vector resonances. 
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1 Introduction 

In the past few years simplified model spectra (SMS) [1-5] have been used as a reasonable 
framework for the interpretation of the limits from mono-jet/mono-photon [6-11] searches 
on direct production of dark matter (DM) at the LHC. In SMS one generally parametrizes 
the cross section in terms of a few parameters, like the couplings of the DM with the visible 
sector, or the mass of the particles assumed to mediate the interaction between the DM and 
the partons in the nucleons, effectively obtaining a way to capture the physical properties 
of a vast class of weakly interacting massive particle (WIMP) scenarios without excessive 
proliferation in the number of free parameters. 

Early LHC results presented by the experimental collaborations were often interpreted 
in the effective field theory (EFT) framework [12-20], which gives a good approximation as 
long as the interaction is mediated by particles with masses well above the collision energy. 
The EFT framework also has the advantage of providing bounds in terms of a common 
contact operator that can be used for comparison with the limits on the spin-independent 
DM-nucleon cross section, obtained in underground experiments like XENONIOO [21] 
and LUX [22], or the spin-dependent cross section, measured for example at PI¬ 

CASSO [23, 24], COUPP [25] and, indirectly, IceCube [26, 27] and ANTARES [28]. 

However, at the center-of-mass energies typically probed in a collider environment it 
is often necessary to consider situations where the approximations underlying the EFT 
framework no longer apply [29-34]. Several recent analyses that compared LHC and direct 
detection (DD) bounds on DM have therefore rather used the SMS framework [35-48]. A 
wide range of possibilities for the DM mediators and couplings has been discussed in [35-48] , 
but in general the literature agrees in that for WIMP masses above ~ 5 GeV the bounds 


- 1 - 





on fjp coming from reinterpretations of mono-jet searches are not competitive with the 
bounds from LUX or XENONIOO. 

On the other hand, because of their reduced number of parameters, the most com¬ 
mon SMS are not intrinsically equipped to capture some of the interesting phenomenology 
of more realistic theoretical DM models (for example those emerging in low scale super- 
symmetry (SUSY)). In models that involve a richer spectrum of particles, several effects 
can arise which are missed in the most simple SMS, like long decay chains, or the well 
known fact that interference between different diagrams, or cancellations in the couplings 
can give rise to “blind spots” for DD searches [49-55]. It is also known that with a larger 
number of particles with differentiated properties one can make a more effective use of the 
complementarity of different experimental strategies, which can be employed in combina¬ 
tion [41, 56-58]. 

In this regard, then, it would perhaps be interesting to give a detailed look at the 
detection issues arising in cases when one moves just one step beyond the SMS approach, 
i.e., when one tries to build models that are halfway in between those SMS characterized 
by just one type of mediator and interaction mechanism, and a UV complete model. We 
take this approach in this paper, in which we combine existing SMS in pairs, with the 
goal to somewhat mimic the behavior of a developed UV theory without at the same 
time drastically increasing the number of parameters, or including the full spectrum of a 
specific model. We only consider SMS and parameter ranges for which the spin-independent 
scattering cross section is substantial, so that a comparison between the limits from DD and 
the LHC is always possible. The combinations involve three popular SMS characterized 
by vector mediators, scalar mediators, and colored scalar mediators, and we take the DM 
particle to be a Dirac fermion. The models given here represent only a few motivated 
examples of the many combinations that can be constructed. We dedicate special attention 
to the blind spots for DD, which stem from interference effects among different diagrams. 
We show that each of the emerging blind spots can be tested in Run 2 of the LHC, by 
different experimental strategies. Our study is thus complementary to previous studies in 
this direction [49-55], although, to the best of our knowledge, the combinations we consider 
here have not been analyzed before in this setting. 

The paper is organized as follows. In Sec. 2 we review the three known SMS of DM 
that constitute the building blocks from which we construct new models for DM searches. 
In Sec. 3 we present the parameter space of these new models and confront them with 
the current bounds from searches for DM, which include underground detectors and DM 
searches at the LHC. We further identify the regions that give rise to suppressions in 
and apply to the resulting blind spots existing bounds from the 8TeV and 13TeV LHC. 
We also calculate the corresponding projected reach of the 14TeV run, showing that this 
can significantly probe these new regions of the parameter space. We finally present our 
summary and concluding remarks in Sec. 4. 
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2 The model blocks 


In this paper we present a phenomenological analysis of a few models of DM that should 
provide “less simplified” model frameworks (LSMS) that to some extent mimic the char¬ 
acteristics of more generic UV theories. 

Model 1. Combining vector and Higgs portal mediators; 

Model 2. Combining Higgs portal and t-channel mediators charged under color; 

Model 3. Combining vector and t-channel mediators charged under color. 

In all three cases we take the DM particle to be a Dirac fermion singlet under the SM 
symmetries. 

The combinations given above are built out of well known SMS, which we call the model 
blocks, whose characteristics we briefly recall in the next subsections. The model blocks 
should respect the gauge symmetries of the SM and not violate minimal flavor violation 
(MFV) [4, 59]. These requirements constrain the allowed forms of the models and their 
parameters, in contrast to the most general forms [2]. 


2.1 Vector mediator 


The mediator considered is a leptophobic Z'. Dark matter SMS based on vector mediators 
of this kind where studied, e.g., in [35, 36, 40, 43, 44, 47, 48, 60-65]. The Dirac fermion 
singlet DM particle, y, is coupled to the new gauge boson, Z'. The new mediator is 
assumed to have negligible mixing with the Z boson of the SM, and to not couple to the 
SM leptons, so that one can easily evade the strong limits from di-lepton resonances at the 
LHC [66, 67]. Note that this assumption will not affect the DD constraints and the results 
of the next sections will remain general. We also always assume that possible anomalies 
are canceled by new heavy states above m,^/, which do not contribute to phenomenology 
at the LHC.^ 

The interaction terms relative to DM detection at the LHC and in underground de¬ 
tectors are 

L D Z'^ja^{g\ - g^-f5)x + Z'^QiYigq - 9fl5)qi , (2.1) 


where the index i runs over the quarks and we have universal vector (axial-vector) quark 
couplings g^ {dq)- The corresponding vector and axial-vector couplings to the DM are g^ 
and g^, respectively. 

In this paper we limit ourselves to the case where WIMPs are produced at the LHC 
through an on-shell mediator: 2my- < mz' ■ In this regime the production cross section 
and mediator width are largely independent of the spin structure of the couplings [5], so 
that we can set either g^^^ or g:^^^ to zero without loss in generality. We limit ourselves 
here to cases of vector boson exchange, which give contribution to the spin-independent 
WIMP-nucleon scattering cross section, <7®^. 

We are thus left with 4 parameters for this simplified model. 


r V v^ 

{m^,mz>,g^,gq I , 

^This assumption is not always warranted, see for instance Ref. [68]. 


( 2 . 2 ) 
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where is the WIMP mass and mz' the Z' mediator mass. This number can be reduced 
to 3 when = 9^, since only the product of the two coupling constants matters for 

2.2 Scalar mediator/Higgs portal 

The second building block is a model such that the fermion DM singlet, Xj is coupled to 
a new singlet real scalar, s. This model has been analyzed, e.g., in [59, 69, 70], but scalar 
mediators have also been studied in, e.g., [45, 46, 71-76]. We note that there are many 
variations with scalar mediators or extended Higgs sectors but this model allows us to keep 
the DM as a SM singlet and allows combinations with the other model blocks. 

The terms in the Lagrangian relevant to DM searches are 

C D-y^XXS - 9ss\<^\‘^ - , (2.3) 


where is the Yukawa coupling between the DM and the singlet, and gs is a mass term that 
induces mixing between s and the SM Higgs doublet, that gives rise to the Higgs boson 
after electroweak (EW) symmetry breaking. In Eq. (2.3) we neglect eventual polynomial 
terms of the singlet only, which are not important for DM SMS (at least at the tree level), 
and which we assume are fixed by the UV completion. We assume that <I> develops the 
SM vacuum expectation value (vev), u: <I> —)• l/\/2(0, u + h)'^ , which can be determined 
in terms of the SM mass and quartic couplings. 

The gs and As Lagrangian terms produce an off-diagonal component in the {h, s) mass 
matrix. The mass matrix is diagonalized by a mixing matrix parametrized by a mixing 
angle 9, 

(hsM] _( cose sme\ I h\ 

\ H J y-sine cosO) \ s J ' ^ ’ 

In general we will identify the lightest scalar with the observed SM-like Higgs. Higgs physics 
measurement, EW precision tests, and vacuum stability then constrain | sin0| [77-80]. 

After diagonalization the relevant terms in the Lagrangian for DM phenomenology are 


C D -y^ {hsM smO + H cos 9) XX 


1 

71 


{hsM cos 9 — H sin 9) 

f 


(2.5) 


where yf are the SM Yukawa couplings and /, / SM fermions. This results in the presence 
of a heavy scalar mediator, H, as well as the SM Higgs, hsM, that couple the DM to the 
quarks. Note that the quarks couple to both mediators in proportion to their Yukawa 
couplings, as a result of the mixing, automatically in agreement with MEV. In the spirit 
of phenomenology one can trade the parameters and gg for 9 and the mass of the heavy 
scalar niH to produce a simplihed model of DM. Note that we implicitly assume that the 
SM and singlet polynomial terms neglected in Eq. (2.3) are such that the ratio of the Higgs 
vevs, tan/3 = v/vg, is much smaller than 1, so to be comfortably inside the bounds from 
perturbative unitarity [80]. 

The DM simplified model is hnally described by 4 parameters. 


{my.,mH,sm29,y^} . 


( 2 . 6 ) 
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Figure 1: Contours of equal spin-independent scattering cross section, in pb in the 
{rriH, Vx) plane for the Higgs portal model. We set 6 = 0.3 and = 100 GeV. The solid 
red line represents the current 90% C.L. upper bound from LUX [22]. The dashed red 
line gives the expected sensitivity of XENON-IT for 2017 or so [81], which we use as an 
example of the generic reach of tonne-scale detectors [82-84]. 


In the next section we will consider, for example, a case where the Higgs portal de¬ 
scribed here interferes with the model of Sec. 2.1. However, we foreshadow this discussion 
by reminding the reader that there is also interference within the Higgs portal model itself, 
as recently pointed out in Ref. [59] . Thus, the contributions to £7®^ intrinsic to this model 
are also going to be affected by interference effects. 

We show this in Fig. 1, where we present contours of in pb in the {mui Vx) plane 
for a fixed 9 = 0.3^ and = 100 GeV. The cross section depends on through the 
reduced mass of the DM proton system and is thus mostly insensitive to changes in 
when 3> m-p. The current upper bound on from LUX [22] is shown as a solid red 
line. It is expected to be improved by future tonne-scale underground detectors [81-84], 
whose reach is shown here as a dashed red line. 

In general the effective interaction between the quarks and the DM will grow with 
the mixing angle, the coupling strength and the mass of the heavy scalar, albeit very 
mildly when mn • When mn ~ ^/ism contributions due to /ism and H cancel 

out and is suppressed. This creates a blind spot for DD, characterized by very little 
sensitivity for underground experiments. 

^ In the analysis of Sec. 3 we will most often assume 6 — 0.2, which is allowed by the bounds from Higgs 
searches at the LHC in the range mn > 125 GeV [80], by perturbativity of the couplings, and one-loop 
corrections to the W mass [80, 85]. As pertains to Fig. 1, one should note that ff — 0.3 is also excluded 
for rriH < 125 GeV by searches at LEP and the LHC [80], and that the figure is here merely intended as 
demonstrative of the effects of interference on the WIMP scattering cross section. 
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2.3 Scalar t-channel mediators 

For the third building block we consider colored scalar mediators coupling the DM directly 
to the quarks. Since these are exchanged in the t-channel for DM production at the LHC 
they are often called t-channel mediators. This case has been analyzed, for example, in [37- 
39, 41]. The new scalars must be charged under color and flavor, so that we can borrow 
the notation used to describe the squarks in the Minimal Supersymmetric Standard Model 
(MSSM), q, even if our model is not necessarily SUSY based. We assume universality 
between the first and second generation for the masses of the “squarks”, as well as their 
couplings to the DM, as required by MFV. For simplicity, for the first two generations 
we further assume that all squark masses and couplings are universal, thus neglecting the 
differences between the left and right components, and up and down squark fields. The 
squarks of the third generation will instead be decoupled and assigned a large mass to avoid 
the complications associated with dedicated stop and sbottom searches. Besides, these will 
not contribute to DD observables. 

The interaction terms then become, 

CD ^ gtj (ulj^xPRUi + u\LXPLUi + d\j^xPRdi + dl^xPidi^ + h.c., (2.7) 

i=l,2 

where Qq is the coupling strength, Ui^L{R) are the left (right) up-type squarks of the ith 
generation, are the left (right) down-type squarks of the ith generation, Ui (di) are 

the up (down) quarks of the ith generation and Pr and Pl are the right and left chiral 
projection operators, respectively. We assume that the stability of the DM is protected by 
a discrete symmetry similar to R-parity. We repeat that all squarks in Eq. (2.7) have the 
same mass in our approximation. 

This simplified model is described by 3 parameters, 

{my^,mq,gq} , ( 2 . 8 ) 

where rriq is the universal squark mass and gq is the universal squark-DM coupling. Note 
that in order to preserve the SM gauge symmetries the squarks will develop the same SM 
gauge interactions as the quarks as happens, e.g., in the MSSM. On the other hand, EW 
gauge and Higgs couplings will not play a role here. This is different from the MSSM, 
where gq is fixed by neutralino mixing and gauge couplings, while in this scenario it is a 
free parameter. 

3 Methodology and analysis of the combined models 

We here confront the LSMS introduced at the beginning of Sec. 2, Model 1, 2, and 3, 
with the bounds from DD of DM in underground detectors and a number of results from 
the LHC: mono-jet searches, searches with jets -|- missing Et (MET), invisible branching 
fraction of the Higgs boson, and bounds on new heavy Z' resonances from the tt and di-jet 
invariant mass distributions. We will in particular highlight different strategies that can be 
used to test the parameter space invisible in DD searches. We also calculate projections for 
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the LHC 14 TeV for each model and compare them to the expected sensitivity of tonne-scale 
underground detectors. 

Models 1, 2, and 3 are implemented in FeynRules [86], which is then used to generate 
CalcHEP [87] and Universal FeynRules Output (UFO) model files. The spin-independent 
cross section is calculated with micrOMEGAs v.4.1.8 [88]. For the mono-jet search we gener¬ 
ate events at the LHC with MadGraph5_aMC@NL0 [89] and PYTHIAS [90] using MLM matching 
up to two jets. For models with squarks we produce both yy+jets, X 9 +jets associated pro¬ 
duction, and gg*-|-jets following the example of [41]. To derive exclusion bounds for the 
mono-jet and jets-pMET searches we use CheckMATE [91-99] and two codes developed by 
some of us, which were previously used in Refs. [100-102]. The LHC 14 TeV projections 
for the mono-jet and jets-pMET searches are obtained by implementing the experimental 
cuts described in [103] and [104] in our own codes and CheckMATE. We do not consider 
other mono-X searches beyond the mono-jet to keep the number of searches considered to 
a minimum and since projections for 14 TeV are not currently available. 

Eor the direct searches for Z' we use MadGraph5_aMC@NLD to compare production cross 
section times branching ratio for Z' —)• tt and Z' —)• qq to the limits given in [105-107]. 
Eor the 14 TeV projections we use [108]. Eor the dijet search we use a combination of the 
8 TeV data for mz' < 1.5 TeV and 13 TeV data for mz' > 1.5 TeV. In the case of the 
13 TeV search the collaboration has recast the search explicitly in terms of the simplified 
model considered here, giving limits in terms of and mz' ■ The collaborations have not 
provided a projection for 300 fb“^ at 14 TeV in terms of Z' —)■ qq, so that we use the di-top 
search throughout the paper for consistency. Einally, we calculate the partial width of the 
SM Higgs to DM particles, with CalcHEP and compare to the limit from [109] 

following the method described in [57]. 

We do not apply to our models the bounds from the relic density of DM. In a sense, 
this is equivalent to considering the most general phenomenological case, neglecting all 
possible bias from the particular thermal history of the early Universe. The reader should 
bear in mind that from this point of view the models analyzed in this section are incom¬ 
plete. In order to satisfy the constraint from the relic density one must consider one or 
more of the following possibilities: large or non-perturbative couplings; additional parti¬ 
cles/interactions; or a non-standard thermal history. 

3.1 Model 1: combining Z' and Higgs portal 

The first LSMS we consider is motivated by the observation that many UV complete 
models with a Z' also contain an extended scalar sector, see for example [HO]. Inspired 
by constructions of this kind, we consider a Z' vector boson associated to a new symmetry 
U (l)x , and a hypothetical extended scalar sector that will include, among others, a U (l)x- 
neutral SM singlet field s that couples to the SM Higgs and the DM particle like in Eq. (2.3). 
If all other degrees of freedom are decoupled, the low energy Lagrangian is just the sum of 
Eq. (2.1) and Eq. (2.5).^ 

different combination of the same two mediators is given, e.g., in [111]. 
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As explained above, we consider only the vector couplings to the Z'. If we make the 
assumption that we are left with 6 free parameters, 


{ 


m 


X’ 





(3.1) 


We assume mixing as maximal as is allowed by the LHC constraints, perturbativity of the 
couplings, and EW precision observables: 9 = 0.2 (see Footnote 2). The cross section 
depends mildly on the angle, via sin 20. 

For a nuclear element N, the differential WIMP-nucleus scattering cross section in the 
non-relativistic limit is given by [112] 


ddyN 

^qF 


-^[Zfp + iA-Z)UfF\Q), 

TTV^ 

X. 


(3.2) 


where jqj is the transferred momentum, Z is the atomic number, A the atomic weight, 
is the average speed of the DM in the halo, and F{Q) is the Wood-Saxon function given 
in [112] as a function oi Q = |q|^/2m7v . 

The contributions to fp and fn come from the effective interactions of the WIMP with 
protons and neutrons, respectively, and they are approximately given in Model 1 by 


fn 


fp 


0xS™20 f mp 

V ml J V 


fTq+‘^ fTG 


, 3 gX 




(3.3) 


where mp is the nucleon mass, and fTq and fxG are the hadronic matrix elements defined 
for example in [112]. We adopt in our study their default values embedded in micrOMEGAs: 
frd = 0.0191, fTu = 0.0153, fxs = 0.0447, and /tg = 1 — fru — frd — frs- 

If > 0, mn m ^^^^, and = <7^ destructive interference between the terms in 
Eq. (3.3) does not take place. In Fig. 2(a) we show contours of in pb in the 9x/q^ 
plane in this case, for fixed values of the DM mass, m.^^ = 10 GeV, and of the mediators’ 
mass, mu = 600 GeV and mz' = 1000 GeV. The mediators’ mass has been chosen so to 
be within the present limits from the LHG for reasonable choices of the couplings, not too 
far from their SM strengths, and also so that we are reasonably far away from interference 
effects between the light and heavy Higgs bosons. 

As was the case in Fig. 1, the solid red line shows the 90% C.L. upper bound on 
from LUX and the corresponding red dashed line the projected reach of XENON-IT. The 
upper bound from mono-jet searches at the LHG 8TeV [9] is shown as a solid purple line, 
and the projected reach for mono-jet with 300 fb~^ at 14TeV [103] is shown with a dashed 
purple line. The solid orange line shows the upper bound on heavy resonances from the 
invariant ti mass distribution after the LHC 8TeV run [105] and the solid cyan line gives 
the equivalent limit in the qq channel [106]. The projected reach for the tt search at the 
LHC 14 TeV [108] with 300 fb“^ does not improve on the 8 TeV bound for mz' = 1000 GeV 
and is consequently not shown here. Finally, the green solid vertical line gives the upper 
limit on the yy_ coupling from a CMS/ATLAS combined analysis of the invisible decay of 
the Higgs boson [109]. 
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Figure 2: (a) Spin-independent scattering cross section in pb in the {y^, 9^/q) plane for a 
model combining Z' mediator and Higgs portal. We set = 10 GeV, mz' = 1000 GeV, 

9 = 0.2 and mn = 600 GeV. The solid red line represents the current 90% G.L. upper 
bound from LUX and the dashed red line is the projected limit for tonne-scale detectors 
in 2017 or so. Solid purple line is the 95% G.L. upper bound from the 8TeV mono-jet 
search at ATLAS [9]. Dashed purple line is the projected limit from the mono-jet search 
at 14TeV with 300 fb~^ [103]. Solid orange line gives the upper limit from searches for 
heavy vector resonances in the di-top channel at 8 TeV [105], and the solid cyan line gives 
the equivalent limit in the qq search [106]. Green solid line gives the upper limit from the 
invisible width of the Higgs boson in a CMS/ATLAS combined analysis [109]. (b) Same 
as (a) but = 100 GeV. (c) Same as (a) but the sign of is negative, (d) Same as (c) 
but ruy- = 100 GeV. 
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As is well known, the DD detection bound on from LUX is significantly more 
constraining then any of the collider limits, even for a relatively light DM mass. LUX loses 
sensitivity with respect to the collider bounds when < 5 GeV. For ra-^ < 62 GeV, the 
invisible width of the Higgs boson analysis places an upper bound on that is stronger 
than the projected reach of many tonne-scale detectors, as was also recently pointed out 
in [57], 

When the WIMP mass becomes larger than 10 GeV the DD bounds become more 
severe, and they reach their maximal sensitivity when ~ 50 GeV. We show the bounds 
for = 100 GeV in Fig. 2(b). Note that the invisible width bound does not apply to this 
case, so that the strongest limits on both couplings here come from DD experiments. Also, 
the projected upper bound on from mono-jet searches at the LHG 14TeV is stronger 
when < 62 GeV, because the XX pair can be produced via an on-shell SM Higgs. 

On the other hand, if < 0, or if it is positive but ^ = — 5 ^, the diagrams corre¬ 
sponding to the Z' and Higgs portal interfere destructively and <7®^ becomes suppressed, 
as can be inferred from Eq. (3.3). Note, incidentally, that the cancellation can only hap¬ 
pen in the nonrelativistic limit. To see this one can consider the relativistic WIMP-quark 
scattering, q{pi)xiP 3 ) q{P 2 )x{Pi)i for which the squared amplitude reads: 

2 _ sin^ 29 ygvliml + pip2){m\ -k p^p^) {IQml - 8pip2)(16m^ - Sp^Pi) 

[{Pi - P2Y - [{pi - P2Y - m?z,Y 

16 sm2e yqy^g^g^ mpm^{pi + P 2 Y {P3 + Pa) t, , . 

\/2 [{Pi - P2Y - ml^J[{pi - P2Y - m%,] 

Obviously, without a detailed knowledge of the incoming and outgoing 4-momenta, a can¬ 
cellation between the terms of Eq. (3.4) cannot take place. Additionally, in Eq. (3.3) there 
is also an effective term due to the coupling of the Higgs to the gluon field strength, which 
determines the position of the blind spot. 

The effects of the cancellation are shown for = 10 GeV in Fig. 2(c). The blind spot 
is in the plots a narrow diagonal region, over which the value of 0 "®^ visibly drops below 
the potential reach of tonne-scale detectors. The corresponding case for = 100 GeV is 
shown in Fig. 2(d). 

From Eq. (3.3) one can derive the condition for the blind spot: 


Vx 


8.22 X 10^ GeV^ 


rrir. 


sin 29 


V V 

9x9q 



(3.5) 


Condition (3.5) shows that the contributions to the amplitude of the diagrams from the 
Z' and Higgs portal are of comparable size for comparable coupling strengths if mz' is at 
least of the order of a TeV or larger. Here we neglect further loop corrections which may 
reintroduce a suppressed coupling to the nucleus. When Eq. (3.5) is satisfied the model is 
beyond the reach of DD searches in underground detectors but can be tested by collider 
means. 

We therefore show in what follows how the interplay of mono-jet searches, searches for 
Z' resonances, and Higgs width measurements introduced above constrains the remaining 
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Blind spot 1: H and Z' 
m)( = \0 GeV 
0=0.2 mH=600 GeV 
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Figure 3: (a) Interplay of LHC constraints for the blind spot parametrized by Eq. (3.5) 
in the (m^/, 9^/g) plane. Here = 10 GeV, mn = 600 GeV, and 9 = 0.2. Solid orange 
line gives the 95% G.L. upper bound from the ATLAS search for heavy resonances in the 
di-top channel at 8TeV. The dashed orange line gives the corresponding reach of the LHC 
14TeV with 300 fb~^ [108]. Solid cyan line gives a combination of the upper bounds on 
the Z' invariant mass in qq searches at 8TeV [106] and 13TeV [107]. Green solid line gives 
the combined ATLAS/GMS upper bound from the invisible width of the Higgs. Dashed 
purple line gives the projected mono-jet reach at ATLAS with 14 TeV and 300 fb~^. (b) The 
bounds projected to the {mz', my^) plane for = 0.2 and rriH = 600 GeV. Gorresponding 
bounds are now lower limits on the masses. The combination of 8 and 13 TeV data in qq 
searches excludes the parameter space within cyan bars. 


parameters of Model 1, when the condition (3.5) for a blind spot is satisfied. In Fig. 3(a) 
the bounds are shown in the {mz', 9^/g) plane, assuming Eq. (3.5) holds. The DM mass is 
here fixed, = 10 GeV, and we have set mn = 600 GeV. Note that the plots essentially 
do not change much for any mn ^ 200 GeV. 

The shaded region at the top of Fig. 3(a) is not allowed, as becomes there nonper- 
turbative, > dvr. In both panels the color code is the same as in Fig. 2. Note that the 
strongest upper bound is currently given by the invisible width of the Higgs, but one must 
remember that the bound does not apply for rriy. > 62 GeV. The limits from the mono-jet 
and direct Z' resonance searches, on the other hand, barely change position over a large 
range of DM masses. The current ATLAS mono-jet bound is too weak and does not appear 
in Fig. 3(a), in agreement with the results of Figs. 2(c) and 2(d). The upper limit from 
Z' —>■ tt is shown with a solid orange line and a combination of the upper bounds on the 
Z' invariant mass in qq searches at 8 TeV [106] and 13 TeV [107] is shown with a cyan solid 
line. Note that the 13 TeV data improved significantly on the 8 TeV data for Z' masses 
above 1500 GeV. 
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Figure 4: A comparison of the limits from mono-jet searches at the LHC and searches for 
Z' in the di-top and di-jet invariant mass distributions in the {g^, g^) plane for mz' = 
1000 GeV, niy- = 10 GeV. 

We calculate that in Run 2 with 300 fb~^ the bounds from resonance searches will 
improve considerably for large masses, as shown by the orange dashed line. When mz' ^ 
1500 GeV the improvement over the 8TeV bound is small and for mz' ~ 1000 GeV we 
project that there will be no improvement in the it resonance search, which explains why 
the orange dashed line does not appear in Fig. 2. Note that the 14TeV mono-jet reach 
is comparable to the current di-top bound at 8TeV and far below the 14TeV projection 
making mono-jet searches less efficient than other strategies for testing this model. 

The mono-jet bounds derived above feature just a mild dependence on the DM mass, 
since the probability for emission of a hard ISR jet is largely independent of it, as long 
as the mass is small compared to the cut on the transverse momentum of the leading jet. 
In Fig. 3(b) we show the parameter space that can by probed by mono-jet searches in 
the {mz', plane for g^/q = 0.2. It can be seen that the mono-jet search is always 
weaker than the direct search for a Z’ in the di-top and di-jet channels. The 14TeV mono¬ 
jet search loses all sensitivity once the DM mass becomes of the order of ~ 350 GeV, as 
expected for the couplings and selection requirements considered here. 

Finally, let us conclude this subsection with a few more comments. Given particular 
values for the Z' mediator mass, the limits shown in the plots are derived under the 
assumption g^ = g^ ■ This choice implies that for most of the parameter space the di-top 
and di-jet searches for heavy resonances are more sensitive to g^jq than the mono-jet search 
for DM. In Fig. 4 we show how the bounds depend on g^ and g^ separately, for values on 
the mediator and DM masses set at mz' = 1000 GeV, = 10 GeV. The case g^ = g^ 
corresponds to the dashed blue line. One can see that Z' searches are more constraining 
than the mono-jet search for a wide choice of couplings, but the limits switch position 
when g^ > 3.5 g^ due to the increased branching fraction Br(Z' —>• xx)- Therefore, one 
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must always keep in mind that even if the bounds presented here are valid over large 
regions of the parameter space, there remains a significant dependence on the underlying 
assumptions. 

For g'^ / the upper bounds in Fig. 2 will move in the (y^, ^9^9x) plane. In 
particular for g^ < g^ the upper bounds from di-top and di-jet searches will become 
stronger, constraining smaller values of the product, and the monojet upper bounds will 
become weaker. For g^ > g^ it will be the other way around. This can be understood 
from Fig. 4: if one takes g^ as large as possible, i.e. ~ 1.4, then the upper bound on 
derived from the monojet search is ~ 0.15, as the figure shows. One can then evaluate 
\/l.4 • 0.15 = 0.45 to get un upper bound on the product. This is essentially the strongest 
the monojet bound can get at the expense of the bounds from the di-jet distribution of the 
Z'. The product is the important quantity for the blind spot and direct detection. 

3.2 Model 2: combining Higgs portal and squarks 

This LSMS features some of the characteristics of SUSY models of DM. In particular the 
bounds can resemble those obtained in cases where the neutralino couples to the SM Higgs 
and additional heavy Higgs bosons; see, e.g., [54]. An obvious difference with the MSSM 
is that in the case presented here the DM is a Dirac fermion with free couplings, and the 
additional scalar is a SM singlet. 

The low-energy Lagrangian of Model 2 is given by the sum of Eq. (2.5) and Eq. (2.7). 
There are 6 free parameters. 


{m^,mq,mH,0,y^,gq]. (3.6) 

As before, we fix 0 = 0.2 because of the mild dependence on sin 20. We start by fixing 
the masses of the mediators, mq = 1000 GeV and mn = 600 GeV. Recall that the latter 
choice implies that we do not investigate interference effects between the heavy and light 
Higgs bosons. 

The effective amplitudes for DD in this case read 


fn 


fp 


sin 28 
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q=u^d^s 


4 


m? — m| 


(3.7) 


where we parametrize the numerical coefficients relative to the tree-level and box-diagram 
WIMP-quark-squark interactions in the proton with Ctree and Cboxj respectively. The 
numerical value of the sum in parenthesis in the second line of Eq. (3.7) is 0.019. Note also 
that we use the micrOMEGAs default constituent value niq = 0.05 GeV for the 3 lightest 
quarks. 

When > 0 cancellations in the amplitude for do not occur. We thus here limit 
ourselves to describing the case y^ < 0, so that interference effects create a blind spot for 
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(a) (b) 

Figure 5: (a) Spin-independent scattering cross section in pb in the {y^, Qq) plane for 
a model combining squark-like mediators and Higgs portal. We set = lOGeV, rriq = 
1000 GeV, 6 = 0.2, and mn = 600 GeV. The solid red line represents the current 90% C.L. 
upper bound from LUX and the dashed red line is the projected limit for tonne-scale 
detectors in 2017 or so. Solid purple line is the 95% C.L. upper bound from the 8TeV 
mono-jet search at ATLAS and solid blue line gives the 95% C.L. upper exclusion bound 
from direct squark searches with jets and missing energy [113]. Green solid line gives the 
upper limit from the invisible width of the Higgs boson from the combined CMS/ATLAS 
analysis. The full parameter space shown is within reach of 14 TeV jets-|-MET and mono-jet 
searches, (b) Same as (a) but m-^ = 100 GeV. 


DD searches. The condition for the blind spot is given by 


Vx 


f 2.05 X 10^ GeV2\ 

V ) 


sin 20 



(3.8) 


We present in Fig. 5(a) contours of in pb in the (y^^, gq) plane for the case = 
10 GeV. The color code for the bounds is the same as in Fig. 2. Additionally, the solid 
blue line shows the upper bound from the ATLAS 8 TeV squark search in jets -|- missing 
Ex [113] (see also [114] for the CMS bound). As was the case for Model 1, the bound on 
\y^\ from the invisible width of the Higgs is stronger than the LUX bound from DD for 
this benchmark point. 

In Fig. 5(b) we show the case with my_ = 100 GeV, for which the bound from the 
invisible width obviously does not apply. As was the case in Fig. 2, the bound from mono¬ 
jet is less constraining for jy^j in this case, as the SM Higgs is not produced on-shell. Thus, 
the parameter space that is not in reach of underground DD experiments remains essentially 
unconstrained, and this is in particular true for the blind spot, when > 62 GeV. 


- 14 - 



















































Figure 6: Solid blue line shows the cross section for squark production through t-channel 
DM exchange at the LHC 8TeV for rriq = 1000 TeV and a WIMP mass = lOGeV, as 
a function of the gq coupling. Solid red line shows the same production cross section when 
= 500 GeV. The equivalent cross sections at the LHC 14 TeV are shown with dashed 
blue and red line, respectively. Solid black line shows the cross section for strong squark 
production at the LHC 8 TeV and the dashed black line the equivalent strong production 
cross section at the LHC 14 TeV [115]. 


On the other hand, we stress that for the squark mass considered in Fig. 5, niq = 
1000 GeV, the 14 TeV jets+MET search with 300 fb~^ at ATLAS [104], as well as the 
14 TeV mono-jet search, are expected to exclude the full parameter space shown in the 
figures. 

To show this in detail, we first present in Fig. 6 plots of the squark production cross 
sections at the LHC as a function of the gq coupling for rriq = 1000 TeV and different 
choices of the WIMP mass. For rriq = 1000 TeV, the cross sections for strong [115] and 
t-channel DM exchange production of the squarks become of equal size when g^ ~ 0.9, and 
this is true at the LHC 8 and 14 TeV alike. 

Then, in Fig. 7(a) we show the dependence of the bounds on the mediator and DM 
masses when Eq. (3.8) holds. The first feature to note is that, as was the case in the 
previous model, the invisible branching fraction of the Higgs yields the greatest constraint 
when < 62.5 GeV, with the exception of the region characterized by small g^ (how 
small it should be to evade the bound depends on rriq) . It is also clear that at 8 TeV the 
jets-l-MET and mono-jet searches dominate in different regions of the parameter space: 
when gq is small direct production of the DM is negligible with respect to the strong 
production of the squarks, whose cross section is shown as a black solid line in Fig. 6. As a 
consequence the jets-|-MET search has the greatest sensitivity and excludes a squark mass 
rriq < 900 GeV. As gq increases direct WIMP production becomes more important and the 
mono-jet bound overtakes the jets-|-MET bound for g^ > 0.3. 

At 14 TeV the reach of the mono-jet and jets-|-MET searches both increase consider¬ 
ably. Taking the current projections of the backgrounds and cuts [103, 104] at face value, 
it appears that the jets-|-MET search will play a more significant role in constraining the 
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Figure 7: (a) Interplay of LHC constraints for the blind spot parametrized by Eq. (3.8) 
in the (m^, Qq) plane. Here = 10 GeV, mu = 600 GeV, and 9 = 0.2. Solid purple line 
gives the 95% C.L. upper bound from mono-jet at ATLAS, and solid blue line the 95% C.L. 
upper bound from jets -|- missing Et- The corresponding 14TeV projections [103, 104] with 
300 fb“^ are given by a dashed purple line and a dashed blue line, respectively. The green 
solid line gives the upper bound from the invisible width of the Higgs in a CMS/ATLAS 
combined analysis, (b) The bounds projected to the {rriq, m^) plane for gq = 0.4 and 
niH = 600 GeV. Corresponding bounds are now lower limits on the masses. 


parameter space of Model 2 than the mono-jet search for Run 2. 

In Fig. 7(b) we show the limits in the (m^, m^) plane. Both the mono-jet and 
jets-l-MET limits show the well known characteristic shapes of squark searches; however, 
the region where the squark and WIMP masses are close to degenerate is not present here, 
as it corresponds to non-perturbative values of the coupling . Thus, the limits on squark 
masses become more robust in the blind spot than in the individual SMS. Also visible is the 
apparent weakening of the limit on from Br(/isM xx) the squark mass increases, 
because decreases as a function of niq. 

3.3 Model 3: combining Z' and squarks 

In this subsection we consider an LSMS designed to mimic a UV completion characterized 
by an additional U{l)x symmetry that remains unbroken down to collider energies (see, 
e.g., [116]). One could imagine that U{l)x is spontaneously broken by the vev of a scalar 
field 4/ charged under U{l)x, giving rise to a light Z' boson which becomes part of the low- 
energy spectrum. The low-energy spectrum also contains scalar particles charged under 
SU{3) whose origin could be SUSY or something else. The details of the UV completion 
are not important for the phenomenology of the LSMS. 
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Constructing a gauge invariant combination of the two existing SMS requires some care 
since there is more than one way Model 3 could be constructed without breaking gauge 
invariance of the full Lagrangian. A first simple way is to imagine that the DM particle 
and the quarks have the same charges under U{l)x, and the scalar colored particles are 
instead U{l)x neutral. In this case the Lagrangian is a straight combination of Eqs. (2.1) 
and (2.7), with a coupling = 9x ~ 9q ^ such that the extra scalars do not couple to 

the Z'. The resulting phenomenological model is described by 5 free parameters; niq, 
mz', 9^/q^ and gq . However, this model cannot develop destructive interference between 
the diagrams with squark exchange and those with a Z' mediator, so that for the purposes 
of this study it is not very interesting. 

Another way of constructing a gauge invariant LSMS out of a combination of the SMS 
with squarks and a Z' is the following, which allows the squarks to have the same coupling 
to the Z' as the quarks, and could be seen as an approximation of a full UV theory involving 
an extended gauge symmetry and a supersymmetric sector. One needs two fermion SM 
singlet DM candidates, ^ and Ci such that ^ is coupled to the Z' like in Eq. (2.1) and C 
is coupled to the squarks like in Eq. (2.7). The symmetry is conserved if the fields are 
charged under U{l)x according to the following table. 





c 

Qi 

%L/R 

U{l)x charge 

+1 

+1 

0 

+ 1 

+ 1 


where we have normalized all charges to 1 for simplicity, and we leave some freedom in 
choice of the coupling constants. As was the case before we assume that any anomalies are 
cancelled by additional states which do not effect the phenomenology described here. The 
low energy Lagrangian can contain the additional terms. 


£ D + h.c., 


(3.10) 


where T is the field that breaks U{l)x when it gets a vev T vq, + ip, with a decou¬ 
pled physical scalar. After the symmetry is broken, ^ and C mix giving rise to two mass 
eigenstates: xi and X 2 which, if we assume <C yiv^, are almost mass degenerate 

with a mass rriy^ = yiVij, and maximal mixing. 

Despite being apparently rather involved, the phenomenological LSMS is characterized 
by only 6 free parameters, 

{m^,mq,mz',g^,gq ,gq} . (3.11) 

To reduce this number one can make the additional assumption g^ = ±g^ = 

Destructive interference between different diagrams arises when g^ = —g^, as 
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(3.12) 


see also Eqs. (3.3) and (3.7). The condition for the blind spot is thus 



(3.13) 
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(a) (b) 

Figure 8: (a) Spin-independent scattering cross section in pb in the {g^ = —g^, gq) 
plane for a combined Z' + squark mediator simplified model. The masses are fixed at 

= lOGeV, mz' = 1000 GeV, and rriq = 1000 GeV. The solid red line shows the 
current upper limit from LUX and the dashed red line the projected limit for tonne-scale 
underground detectors. Solid purple line is the upper limit from the 8 TeV mono-jet search 
at ATLAS and solid blue line is the upper limit from the 8 TeV jets-|-MET search at 
ATLAS. Solid orange line is the ATLAS limit on Z' resonances from the invariant ti mass 
distribution at 8 TeV, while the solid cyan line gives the equivalent bound from the qq 
distribution. The full parameter space shown is within reach of 14 TeV jets-|-MET and 
mono-jet searches, (b) Same as (a) but = 100 GeV. 

We plot in Eig. 8(a) contours of in pb in the 9q) plane for = 10 GeV. 

The masses of the mediators are here set at rriq = 1000 GeV and mz' = 1000 GeV. The 
color code is the same as in the previous figures. The case with = 100 GeV is shown 
in Fig. 8(b). One can see again that the collider bounds barely move by changing the DM 
mass, but DD bounds reach their close-to-maximal strength when rriy. = 100 GeV. 

Note that the mono-jet and jets-|-MET ATLAS searches yield very comparable bounds 
for this choice of mediator masses. Note also that, for these mediator masses, the 14 TeV 
projected reach with 300 fb~^ for both searches covers the full parameter space, so that the 
limit does not appear in the figure. 

The dependence of the collider bounds on the mediators’ mass when we confine our¬ 
selves to the blind spot for DD experiments (i.e. Eq. (3.13) holds) is shown in Fig. 9. In 
Fig. 9(a) we show the interplay of the collider bounds on the {mz', rriq) plane for g^^^ = 0.1. 
The figure shows that a squark mass below ~ 900 GeV is excluded by the jets -|- missing 
energy ATLAS search, independently of the couplings, because squarks are produced via 
strong interactions. On the other hand a Z' mass below ~ 700 GeV is excluded by the 
di-top search for heavy resonance when g^^^ = 0.1. The line moves to the right for larger 
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Figure 9: (a) Interplay of LHC constraints for the blind spot parametrized by Eq. (3.13) 
in the (m^/, rriq) plane. Here rriy^ = 10GeV, = —g^ = 0.1. Solid purple line gives 
the 95% C.L. lower bound from mono-jet at ATLAS, and solid blue line the 95% C.L. 
lower bound from jets -|- missing Et- The corresponding 14TeV projections with 300 fb~^ 
are given by a dashed purple line and a dashed blue line, respectively. The solid orange 
line gives the current 95% C.L. lower bound on heavy resonances from the di-top distri¬ 
bution. The corresponding projected reach at LHC 14TeV with 300 fb~^ is shown with a 
dashed orange line. The solid cyan line gives the bound on Z' from the qq invariant mass 
distribution at the LHC 13TeV, 4fb“^. (b) Same as (a) but g^ = —g^ = 0.4. 


couplings, as is shown in Fig. 9(b), where we present the case with g^^^ = 0.4. 

The reach of the mono-jet and jets-|-MET searches at the LHC 14TeV, shown as a 
dashed purple and dashed blue line, respectively, are very sensitive to the size of g^^^. 
In general, our projections, based on [103, 104], show a significantly greater reach for 
the jets-l-MET search, which for g^^^ = 0.4 can probe the squark mediator mass up to 
~ 2000 GeV. 

Conversely, the reach of the di-top search for heavy resonances will allow one to probe 
the Z' mass up to more than 3000 GeV for this coupling, showing again, that a comple¬ 
mentary use of different detection strategies can constrain a large part of the parameter 
space invisible in DD experiments. 

4 Summary and conclusions 

In this paper we have considered LSMS as simple extensions to dark matter simplified 
models, with the goal of mimicking some of the characteristics of more realistic models 
without introducing an excessively large number of parameters. Our starting assumption 
is that, to some extent, many of the features of complex models that cannot be explored in a 


- 19 


























single SMS framework, like interference of different diagrams that result in suppressions of 
the DM scattering cross section, or complementary use of different experimental strategies 
to explore a variegated parameter space from different angles, can instead be employed in 
LSMS built from simple combinations of existing SMS. 

We thus considered three cases characterized by a Dirac fermion WIMP coupled to 
more than one mediator: 1) heavy vector mediator and Higgs portal; 2) squark-like medi¬ 
ator and Higgs portal; and 3) squark-like mediator and heavy vector mediator. We have 
limited ourselves to the analysis of SMS characterized by a sizeable spin-independent scat¬ 
tering cross section, and to a range of DM masses for which a direct comparison of 
the limits from the LHC with the limits from DD in underground laboratories is possible, 
i.e., we have not considered cases with < 10 GeV. In each case we have confronted 
the models with the present bounds from DD at LUX, and a number of constraints from 
Run 1 at the LHC: mono-jet searches, direct searches for heavy resonances in the di-top 
and di-jet channels, invisible width of the SM Higgs boson, and searches for colored heavy 
particles in jets with missing Et ■ 

In particular, we dedicated special attention to the parameter space showing interfer¬ 
ence between different diagrams, which gives rise to blind spots for DD experiments. For 
these blind spots we have investigated to what extent interplay between collider limits on 
different particles and couplings can be used for a full exploration of the parameter space 
and we have additionally considered projections for the LHC 14TeV run. 

As is well known, in general for DM mass ranges such that > 10 GeV, the LUX 
upper bound on constrains the coupling constants of WIMP SMS by at least one order 
of magnitude more strongly than any of the LHC searches considered here. Even more 
strikingly, we find that the present bound from LUX also outperforms projected bounds 
for 300 fb“^ at the LHC 14TeV in mono-jet searches and the projected reach of searches 
for heavy vector resonances in most cases. 

A few exceptions to this rule can however be found. One is that in models with a 
Higgs portal, for < l/2m/ig^ the coupling between the Higgs bosons and the WIMP 
are also strongly constrained by the invisible width of the SM Higgs boson, in agreement 
with what was observed in [57]. The other exception, which comprises the central part of 
this work, is that in regions of the parameter space where is suppressed by interference 
effects a combination of LHC searches can effectively place strong limits, especially at the 
end of Run 2. 

In particular, we found that the following features apply to detections of a blind spot: 

• The model involving a Z' and Higgs portal (Model 1) is at present not constrained 
at all by mono-jet searches in the blind spot if = g'^. Moreover, under this 
assumption, the projections for the LHC 14TeV show that searches for heavy Z' 
resonances will constitute the most effective strategy, among the ones considered here, 
to probe this part of the parameter space. Mono-jet searches, even with 300 fb“^ in 
the 14 TeV run, may be competitive with searches for Z' only for cases with g^ ^ g^ . 

• In the two models (Models 2 and 3) involving squark-like mediators the bounds 
from mono-jet and jets-|-MET searches on the coupling gq are at present comparable. 
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However, according to the collaborations’ simulations [103, 104], the constraining 
power of jets+MET searches at the 14TeV LHC for the blind spots significantly 
outperforms the expectations for mono-jet searches. 

In general we find that the complementarity of different search strategies is crucial 
to obtaining the best constraints on DM at the LHC. Many motivated models are not 
necessarily constrained by direct detection experiments and thus demand careful attention 
at colliders. In this paper we have shown some first developments towards describing such 
scenarios in terms of a less-simplified model framework that can be investigated at the 
LHC. 
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